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Abstract 

We present a coalescent process where three particles merge at each coagulation 
step. Using a random walk representation, we prove duality with a fragmentation 
process, whose fragmentation law we specify explicitly. Furthermore, we give a sec- 
ond construction of the coalescence in terms of random binary forests and study 
asymptotic properties. Starting from particles of unit mass, we obtain under an 
appropriate rescaling when tends to infinity a well-known binary coalescence, 
the so-called standard additive coalescent. 

Subject classifications: 60J25; 60J65. 
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Introduction 

Generally speaking, a stochastic coalescence is a Markov process describing the coagu- 
lation of particles characterized by their size only. The rate at which particles merge 
depends just on the members involved. Conversely, fragmentation processes describe a 
Markovian evolution of particles which split independently into new particles (branching 
property). The goal of this paper is to study the stochastic coalescence with ternary 
coagulation kernel 

K(r, s,t) = r + s + t + 3, r,s,t >0, 

to which we will simply refer to as ternary coalescence or ternary coalescent process. 
Here, three particles of sizes r, s, t coagulate into a new particle of size r + s + t at 
rate r + s + t + 3. Although at first glance, the kernel k may look somewhat arbitrary, 
the corresponding process enjoys rather interesting properties. Similar to the additive 
coalescence, that is the coalescence where two particles with masses s, t merge at rate 
k(s, t) = s + t, the state chain of the ternary coalescence admits different representa- 
tions. In the spirit of Bertoin we show how it can be obtained by looking at excursion 
intervals of a one-dimensional conditioned random walk. As a by-product of our repre- 
sentation, we establish duality with a fragmentation process via time-reversal. We stress 
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that this is a unusual feature, because the branching property normally fails when time 
is reversed in a coalescent process. Section 7 of Bertoin gives a brief overview over 
cases where such a duality relation has been proven. See also Chapter 5.5 in Pitman's 
lecture notes [19] for further discussions. 

Using the same construction, we study asymptotic properties of the ternary coales- 
cence starting from particles of unit mass. Properly rescaled in space and time, we 
observe in the limit N ^ oo the so-called standard additive coalescent, which has been 
obtained by Evans and Pitman in [12j as the limit n — t- cxo of the additive coalescence, 
started at time —(1/2) Inn with n atoms of size 1/n. Here, Bertoin's characterization [1] 
of the dual fragmentation process connected to the standard additive coalescent by time- 
reversal plays a pivotal role. We emphasize that even though k is a ternary coagulation 
kernel, we end up in the limit with a binary coagulation process. 

We also highlight a second construction of the ternary coalescence involving random 
binary forests, following the ideas of Pitman in [ISJ . We point out that this representation 
could instead be used to work out our results. 

The rest of this paper is organized as follows. In the first section, we describe the 
semigroup of the ternary coalescence and derive some further properties. We finish 
this part by computing the one-dimensional statistics for the underlying state chain 
starting from an odd number of particles of unit mass. Its special form already hints 
at a connection to hitting times of a one-dimensional nearest neighbor random walk, 
which we elaborate in the next section. There we prove duality via time-reversal with 
a fragmentation process, using an explicit construction of the coalescence in terms of 
ladder epochs. In the third part, we turn our attention to random binary trees and find a 
second representation of the ternary coalescence which is based on random binary forests. 
Finally we use again the random walk representation to study asymptotic properties of 
the coalescence in the last section. 

1 Some basic properties 

Throughout this text, let 

N = {1, 2, . . .}, Z = {. . . , -1, 0, 1, . . .}, Z+ = N U {0}. 

The coalescent process will take values in the space of decreasing numerical sequences 
with finitely many non-zero terms 

= {s = (si, S2, ■ ■ ■) : Si > S2 > ■ ■ ■ > 0, Sk = for k sufficiently large} . 

We may think of elements of a sequence s G iS^'' as (sizes of) atoms or particles and 
simply identify s with its non-zero components. If we write s = (si, . . . ,S;), the non- 
zero components of s are precisely given by Si,...,S/. If s = (si,S2,...) € and 
1 < i < j < k, we use the notation s*®-'®'^ for the sequence in obtained from s by 
merging its ith, jth and kth terms, that is one removes Si, Sj, Sk and rearranges the 
remaining elements together with the sum Si + Sj + Sk in decreasing order. 

Let us define the object of our interest. Recall the kernel k, from the introduction. 
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Definition 1.1. The ternary coalescence with values in S'^ and kernel k is a continuous 
time Markov process X — {X{t), t > 0) with state space Sj' for an appropriate subset 
sj' of S-^, and jump rates 

g(s, •)= ^ K,{si,Sj,Sk)Ssi@j@k. 

l<i<j<k,Sk>0 

This definition can be adapted in an obvious way to other coagulation kernels, leading 
to different stochastic coalescence models, for example the additive coalescence with 
kernel k{s, t) = s + t. 

Before looking at concrete realizations, we collect in this section some basic properties 
which can be read off from the kernel k. and the very definition of jump-hold processes 
of the above type. Denote hy X = {X{t), t > 0) the ternary coalescence, started from 
a finite configuration r = (ri, . . . ,r7v) € S^, where N = 2n -\- 1, n E Z_|_. We write 
M = ri + . . . + for the total mass in the system. For every k = 0, . . . , n + 1, let 
Tfe be the instant of the kth coagulation, with the convention Tq = 0, T„+i = oo. The 
state chain or skeleton chain X' of the coalescent process is given by XI. — X(Tk), 
k — 0, . . . ,n. We use the expression for the number of particles at time t, whereas 
J{t) = max{A; G : < t} stands for the number of jumps up to time t. Note that 
#(t) = iV-2J(t). 

1.1 State chain and semigroup 

Proposition 1.1. In the preceding notation, the following holds true. 

(i) The sequence Ak — Tk — T^-i, k = l,...,n, of the waiting times between two 
coagulations is a sequence of independent exponential variables with respective pa- 
rameters 

a{k) = 1-{M + N + 2-2k){N + l-2k){N - 2k). 
In particular, the sequences {T'fc}o<fe<n o-nd {'^fc}o<fe<n o-f^ independent. 

(ii) The sequence {A'^}o<fe<n 'is a Markov chain with transition probabilities 

Si + s,- + Sfc + 3 



P [Xl^^ = 8^®^®'= I XI = s) 



a{l + l) 



where 0<l<n,l<i<j<k< N — 21, and s = (si, . . . , Sn^2i) E is a generic 
finite configuration with total mass Si + . . . + SAr-2/ = M such that P(A'/ = s) > 0. 

Proof: Let < I < n, and put L — N — 21. By construction, the time A^+i between the 
lih and the {I + l)th coagulation given XI — s — {si, . . . , sl) is exponentially distributed 
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with parameter 



{si + Sj + Sfe + 3) 

l<i<j<k<L 

^ L L L ^ 

{si + + Sfc) - 3 ^ Si - 3 ^ (2si + Sj) 



= ^(M + L)(L-l)(L-2) = a(Z + l). 

Therefore, the waiting times {Ak}k=i,...,n do not depend on the states {A'j^}fe=i,...,n. The 
rest follows from the construction of our process. □ 

We turn to a description of the semigroup. Recall that X starts from A'(O) = r = 
(ri,...,riv). 

Proposition 1.2. In the notation above, consider a partition tt o/{1, . . . , N} into N — 21 

(non-empty) blocks Bi, .... i?jv-2/; each of odd cardinality. Denote by A'^{N — 21) the 
event that the N — 21 atoms of XI result from the coagulation of particles {rj : i e Bj}, 
j = 1,...,N -21. Then, with r^. = XlieSj 

,| N-21 

Proof: The first coagulation involves three particles with labels in the block Bj with 
probability 

^ a(l) 2a(l) 

Now consider an arbitrary sequence [ki, . . . ,ki) taking values in {1, . . . ,N — 21} such 
that for every j = 1, . . . , iV - 2/, \{i < I : h = j}\ = {\Bj\ - l)/2. Using the Markov 
property of X', we see that the probability that for all i = 1, . . . , / , the ith coagulation 
affected only particles formed from initial particles with labels in Bk. equals 





5,|)/2)!(|5, 


-1)! 


((r5^. + l)/2)!((|S, 




l)/2)! 



1 (K. + |i?,|)/2)! 



«(l)---«(Oy (K + l)/2)! 
Observe that the number of such sequences {ki, . . . ,ki) is 



S\B,\ - l)/2, . . . , (|5jv-2/| - l)/2j {{\B^\ - l)/2)\ ■ ■ ■ {{\B^-2i\ - l)/2)\ 
This proves the statement. □ 
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In the setting of the proposition, denote by A^^it) the event that X{t) has N — 21 
atoms, each resulting from the merging of {rj : i G Bj}, j = 1, . . . , N — 21. Since the 
sequence of coagulation times and the skeleton chain X' are independent, 

P (A^(t)) = P (T^ < t < T^+i, A;(iV - 2/)) = P = N-2l)F (A;(iV - 2/)) . 

In particular, the semigroup of X is described by the preceding proposition and the 
distribution of the number of particles at time t, which is computed in the following 
lemma. 

Lemma 1.1. In the notation above, for I = 0, . . . ,n and t > 0, 
Proof: We use 

P =N-2l) = ¥ (T^+i >t)-F{Ti>t). 

Note that is distributed according to Yl'i=i where a{i) is as in the statement 

of Proposition II. ![ and 61,62,... is a sequence of independent standard exponential 
variables. As a general fact, a sum of k independent exponential variables with pair- 
wise distinct parameters a{i) > follows the hypoexponential distribution, that is the 
probability distribution with density 

k k 

Integrating the density and regrouping terms result in the statement of the lemma. □ 



1.2 The monodisperse case 

We turn to the situation where X{0) = r = (1, . . . , 1), that is the coalescent process 
is started from the monodisperse configuration consisting of = 2n + 1 atoms of unit 
mass. In this case, the total mass M equals A^, so the rates a{i) simplify to 

a{i) = {N + l-i){N + l-2i){N -2t). (1) 

If s = (si, . . . , Sm) G S'^ we denote by 7(3) the number of different m-tuples that can 
be built from the elements Si (recall that by our convention Sj > 0). To put it into a 
formula, if {si-}i<i<p is a maximal family of pairwise disjoint non-zero elements from the 
sequence s, and ki = \{j = 1, . . . ,m : Sj = Si-}\, we define 



^ki, . . . , kp) 
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In other words, the ranking map 

oo 

rk: y N"^ 5^ 

m=l 

which orders {ri, . . . ,rm) G decreasingly satisfies |rA;~-^(s)| = 7(s) for each s ^ 
(0, . . .) G S^. As a corollary of Proposition II. 2[ the one-dimensional statistics for A" 
look as follows. 

Corollary 1.1. Let < I < n and s = (si, . . . , Sn-2i) G S'^ with s-i G N odd for all i, 

and Si + . . . + Sn-21 = N. Then, in the situation described above, 



N-2l\l J J-J- Si\^ 

Proof: The the starting configuration is given by (ri, . . . , rjv) with rj = 1 for each i. 
Thus, if X[ has N — 21 atoms of the sizes Si > . . . > SN^2h then there is a partition vr 
of {1, . . . , A^} into N — 21 blocks Bi, . . . , Bn^21 of cardinality \Bj\ = sj, such that the 
atoms of A"/ evolved from merging the particles {r^ : i G Bj}. Denote this event by 
A;(iV-2/). Since 

N\{N-1)\ 
"^'^■■■"(^)=(Ar-/)!(Ar-2/-l)!' 
we obtain from Proposition 11.21 (note that here r^^ = \Bj\ = Sj) 

-1 



P(A;(iV-20).<^-'-^')'^^ 



(A^-1)! 

The number of such partitions vr is given by 

7(s) / N 



n(.-)<4.). 



{N -2l)\\s,,..., 

By multiplying the last two expression together, we arrive at the stated expression. □ 

As the reader may already check at this stage, XI has the same distribution as the 
decreasingly ranked sequence oi N — 21 independent copies C,i of the hitting time 
of a simple random walk, conditioned on + . . . + ^Ar_2/ = A^, and filled up with zeros 
(see Section [2l3] for a definition of these quantities). Indeed, if C,{k) denotes the A;th order 
statistic of ^1, ... , ^n-21, then for s = (si, . . . , Sn-2i) G 

P ((e(iV-2/), • • • , ^(1)) = {SU • • • , SN-2l) I ^1 + • • • + ^N-.2l = N) 

= 7(S) P ((6, . . . , ^N-2l) = (Sl, . . . , SN-2l) I 6 + • • • + ^N~2l = N) , 



and an application of Lemma 12.11 affirms that the last expression coincides with the 
one obtained in the corollary. The connection between random walks and the ternary 
coalescence will become much clearer in the next section. 
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2 Duality with fragmentation via random walks 

Our intention of this section is to prove duality of the ternary coalescence with a frag- 
mentation process. Let us begin with an informal description of such processes. 

Conversely to the phenomenon of coagulation of particles, one often observes in 
nature or science processes of fragmentation. In these systems, particles are broken 
into new ones as time passes. As an example, one may think of DNA fragmentation 
in biology or fractures in geophysics. Just as for coalescent processes, one needs to 
impose constraints on such systems to make them mathematically tractable. First, one 
assumes that the process has no memory in the sense that the future does only depend 
on the present state and not on the past. Second, one supposes that a particle is entirely 
characterized by its size, that is by a real number, and third, one requires the system to 
fulfill the branching property, which means that particles split independently. 

Naively, one might first guess that a coalescent process can always be turned into a 
fragmentation process by running the latter backward in time. However, even though 
the memoryless property is preserved under reversing the time, the branching property 
is normally not fulfilled. In fact, there are only few examples known where a duality 
relation holds (see [6j Section 7 for an overview). 

In view of our informal characterization, it is natural to call a Markov process with 
values in a ternary fragmentation process, if each particle splits at a certain rate 
according to some dislocation law into three smaller pieces, where both the rate and the 
dislocation law depend only on the particle size s, and the sizes of the newly formed 
elements sum up to s. Ranked in decreasing order, these three particles together with 
the ones that did not split form the next state of the process. In particular, different 
particles split independently. 

For our ternary coalescence starting from iV = 2n + 1 atoms of unit mass, we shall 
prove 

Theorem 2.1. Reversing the coalescence chain {A'^}o<fc<n in time results in the state 
chain of the fragmentation process, whose dynamics are given in the forthcoming Propo- 
sitionlKM 

We will derive our result from an explicit construction of the skeleton chain X' in 
terms of (lengths of) excursion intervals of a conditioned random walk. This represen- 
tation will also be useful for studying asymptotic properties in the last section. 

2.1 From configurations to paths to mass partitions 

We first show how subsets of {0, 1, . . . , 2n} can be identified with certain paths of nearest 
neighbor walks on Z of length 2n + 1. The excursion intervals above two consecutive 
(new) minima of such paths partition the space Z/(2n + 1)Z into discrete arcs. Taking 
the ranked sequence of their lengths, we obtain the main object of our interest. 

To begin with, define the configuration space C„ to be the set of all subsets of 
{0, . . . , 2n} which have cardinality less or equal n. We often represent x E hy the 



2 DUALITY WITH FRAGMENTATION VIA RANDOM WALKS 



8 



vector {x{i))o<i<2n, where 



x(i] = 



1 , i E X 
, i ^ X 



Under this identification, we may regard x as a mass distribution. We use the terminol- 
ogy that a site i is occupied by a mass if x{i) = 1 and vacant otherwise. The number of 
occupied sites (the cardinality of the subset x) is denoted by 



We identify configurations x G C„ with paths of a nearest neigbor walk of length 2n + 1 
on Z in the following way. Starting from the origin at time zero, the walk goes one step 
above if site is occupied, i.e. x{0) = 1, and down otherwise, then above if x{l) — 1, 
down if x{l) — and so on, up to time 2n. More precisely, if we let ^i{x) — l{x{i) — 1} 
for < i < 2n, the corresponding path S{x) is given hy S{x)o — and for 1 < j < 2n+l, 



Notice that by definition, S{x)2n+i — 2(|x| —n) — 1. Clearly, the mapping Cn ^ x ^ S{x) 
is one-to-one. 

As we show next, the excursion intervals of such a path provide us with an element 
(pi{x) in the space of cyclically ordered partitions of Z/(2n -|- 1)Z into discrete arcs, 

■Psn+i — — (^1' ■ ■ ■ ' ^rn) '■ there exist ai < a2 < ■ ■ ■ < am < 2n + 1, 
m, Oj e N, such that for 1 < i < m — 1, Sj = [oj, Oj+i) fl N, 
= {[am, 2n + 1) U [0, ai)) n Z+} . 

Take x e Cn, and let M — 2{n — \x\) + 1. With m{x) = mino<j<2n+i 'S'j(a;), define the 
first time at which S{x) reaches rn{x) + k, /c = 0, . . . , M — 1, 



For i — 1, . . . , M, put ai — mM-i{S{x)). We construct a sequence s° = (si, . . . , sm) £ 
^2n+i by setting = [a^, 0^+1) nN for i = 1, . . . ,M-1, sm = ([aM,2n + 1) U [0,ai))nZ+, 
and define v'i(.x) = s°. In other words, if we look for /c = 0, . . . , 2n at the shifted path 
9k{S{x)) defined by 



then the element s° corresponds to the M successive excursion intervals of 9rnj^j_^S{x) 
above two consecutive (new) minima. The length |sj| of such an interval is also referred 
to as a ladder epoch. Defining (fi2 as the function which sends s° = (si, . . . , s^) e 7^2n+i 



x\ = 



|{^e{0,...,2n}:x(^) = l}|. 




mk{S{x)) = inf {j > : Sj{x) = m(a;) -|- k} . 
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Figure 1: The black dots represent the configuration x = {0, 4, 7, 8, 9, 11, 12} C Cg- The 
corresponding path S{x) starts at zero and ends in —3 at time 17. It is periodically 
extended up to time 20 to better recognize the excursion intervals (pi{x). They are 
visualized by the dashed line above the x-axis, where the crosses mark the endpoints of 
the intervals, i.e. ipi{x) — (si, S2, S3) with si = [3, 4), S2 = [4, 7), S3 — [7, 17) U [0, 3). 



to its arc lengths {|sj|}i<j<jji, arranged in decreasing order, we finally obtain an element 
in the space of mass partitions 

s = (si, . . .,Sm) : si> S2> . . ■> Sm,m,Si eN,^Si^2n+l 

i=i 

By filling up with an infinite sequence of zeros, we will often identify mass partitions 
with elements in S^. To summarize our construction, the concatenation map ip 

sends configurations a; e C„ via their path representations to partitions of Z/(2n + 1)Z 
and then to mass partitions. 



2.2 Random evolution 

Our purpose here is to randomize the input of the map ip : Cn ^ ^2n+i obtain (a 
sequence of) random mass partitions. More precisely, we construct two Markov chains 
on Cn running from time zero up to n as follows. Let X — {-'^fe}o<fe<n be the Markov 
chain with Xq — and transition probabilities 



X G y and y\x = {i} for some i e {0, . . . , 2n}\x 
otherwise 



In words, {Xq, . . . , Xi) is obtained by occupying successively / sites from {0, . . . , 2n}, 
chosen uniformly at random. From the point of view of sets, Xi is uniformly distributed 
on the space of all Z-subsets of {0, ... , 2n}. By identifying with the random path S{Xi), 
we will also think of Xi as simple random walk up to time 2n + 1, conditioned to end at 
position —2{n — I) — 1. 
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Let Y = {Yk}o<k<n be the Markov chain with Yq being uniformly distributed on the 
space of all n-subsets of {0, . . . , 2?7,} and transition probabilities 



PY{x,y) 



1^ , y G X and x\y = {i} for some z G {0, . . . , 2n}\y 
, otherwise 



In words, {Yo, . . . ,Yi) is obtained by removing successively / masses chosen uniformly at 
random from the starting configuration Yq. In terms of sets, Yj is uniformly distributed 
on the space of all {n — I)- subsets of {0, ... , 2n}. Similarly as above, Yi can be identified 
with simple random walk up to time 2n + 1, conditioned to end at —21 — 1. Note that 
by construction, we have the duality relation 

(Xo, . . . , Xn) = {Yn, . . . , Yo). (2) 



2.3 Realization of the skeleton chains 

We are not interested in X and Y themselves, but rather in ip{X) = {f{Xk)}o<k<n and 
<f{Y). As we will show in Proposition 12. 3^ the former is the state chain of the ternary 
coalescence starting from = 2n + 1 atoms of unit mass. The latter is characterized 
by Proposition 12.21 as the state chain of a fragmentation process starting from a single 
particle of mass A^. 

We need some preparation. Recall that simple random walk on Z is the Markov chain 
S = {'S'm}m>o with So = and 5*^ = Ci + • • • + Cm, where d, (2, ■ ■ ■ are independent 
random variables with W{Q = ±1) = 1/2. For A; G Z, the first hitting time of k is 
denoted by 

Hk = inf{m > 1 : Sm = k}. 
The following result about the distribution of Hk is classical. 

Lemma 2.1. Let A; G Z, k ^ 0, and m G N. Then 

™/rr N f -ii ^?ih/9)2~"' , A; = m[mod2l 

P {Hk =m) = < V(m+|fc|)/2; ' L J _ 

1 } k ^ m[mod 2] 

Moreover, if m = 2n + 1 and k is a fixed odd number, as n ^ 00, 



F{Hk = m)^l 



3' 



2 V vrra 

Proof: Clearly, for the probability to be different from zero the numbers k and m must 
have the same parity. Then, using the hitting time theorem (see for example |14ij) in the 
first equality. 

The second statement follows from Stirling's formula for the factorial. □ 
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Before looking at (p{X) and ^p{Y) in detail, let us give an indication that the former 
is the skeleton chain of the ternary coalescence. Recall Corollary 11.11 and the connection 
between X' and hitting times. Let N = 2n + 1, < I < n, and take independent 
copies C,i of the hitting time Denote by C,{k) the kth order statistic of ^i, . . . , C,n~21- 

Proposition 2.1. (f{Xi) is distributed according to {^{n-21), ■ ■ ■ conditionally on 

^1 + ■ ■ ■+^N-2i = N, i.e. the one- dimensional distributions of (f{X) and X' started from 
N atoms of mass one agree. 

Proof: We identify Xi with simple random walk S{Xi) up to time A^, conditioned to 
end at —[N — 2/). For notational simplicity, let us write S instead of S{Xi). Also recall 
the definitions of ^^('5') and mfc(S') from Section 12.11 By Theorem 1 of [8], if is a 
uniform random variable on {0, . . . , A^ — 2/ — 1} independent of S, then the chain {S) 
has the law of S conditioned on iJ_(jv-2Z) = X. Moreover, the index m,^ is uniformly 
distributed on {0, . . . , A^ — 1} and independent of the chain 6m„ (S). Denote by dkXi the 
shifted configuration defined by 9kXi{i) = Xi{i + /c[mod A^]). Clearly, ^{Xi) = (p{6kXi) 
for each k. From Theorem 1 of [8] we thus infer that for {si, . . . , Sn^2i) 

P {^{Xi) = (Si, . . . , SN-2l)) = P {fiOn.,M = (Sl, • • • , SN^2l)) 

= P ((^(^-20, • • • , ^(1)) = (Sl, . . . , SN^2l) Kl + • • • + ^N-2l = N) . 

□ 

For the moment, we leave (p{X) aside and first turn to ^p{Y). In the sequel it is 
convenient to use the notion of multisets, which we distinguish from normal sets by 
using double braces. For example, {{a,b,c,c}} contains the elements a, b each with 
multiplicity 1 and the element c with multiplicity 2. The cardinality of this multiset is 
4, the order of elements is irrelevant, as for sets. 

Let ^1,^2,^3 be three independent copies of the hitting time H_i. To state the 
transition mechanism of ip{Y) in a concise way, we define a family /i = {fig, s > 3 odd) 
of probability laws, supported on 

= {R= {{ri, r2, rs}} : G N odd, n + r2 + = s} , 

by setting 

fXsiR) = P ({{6,6, 6}} = I 6 + 6 + 6 = . (3) 



More explicitly, applying Lemma 12.11 results in the expression 

-1 



/ D\ ^ ri \ r2 \ ri 

MR) = 7 



ri+l ) \ r2 + l j \ rs+l 
2 / \ 2 / V 2 



s+3 
2 



(4) 



3rir2r3 

where 7 is the number of triplets {ri,rj,rk) that can be formed from R = {{ri,r2,ri}} 

6 , |{ri,r2,r3}| = 3 
7 = "( 3 ' I{'"i,^2,r3}| = 2 . 
1 , \{rur2,r3}\ = 1 
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Proposition 2.2. (p{Y) = {<^(yfc)}o<fc<n is a Markov chain. Its transition mechanism 
from time I < n — 1 to I + 1 is described as follows. 

(a) Conditionally on (p{Yi) = s = (si, . . . , S21+1) G 'P^n+u select an index i e {1, . . . , 21+ 
1} according to the law 

P(. = ,|^(y,) = .) = |^. 

(b) Given ^p{Yi) = s and l = i, split Si according to the law fig, into three numbers and 
rank them together with Sm, m G {1, . . . , 2i + l}\{i}, in decreasing order to obtain 
a new mass partition. 

Proof: Fix / G {0, . . . , - 1}. We write ip{Y)o.,i for tlie vector {(p{Yo), . . . , (p{Yi)). The 
Markov property will follow from 

(i) (p{Y)q.i and (f{Yi+i) are conditionally independent given (pi{Yi). 

(ii) (fii{Yi) and (p{Yi^i) are conditionally independent given (p{Yi). 

Indeed, assuming (i) and (ii), we have for tq-i+i = (fq, . . . , r^+i) e P2n+i x • • • x ^2n+i' 

P((^(r)o:;+i = ro:,+i) 
= J2 ^Hy)o:i = ro, I ipiiYi) = u°)P((^(y,+i) = r,+i | ip,{Yi) = u°) 

u:v32(u°)=r; 

X P((^i(l^) =U°) 

- P ((^(ll+i) = r,+i I ip{Yi) = ri) P ((^(F)o, = ro,) . 

For (i), the key step is to show that the conditional law of (p{Y)q.i given Yi only depends 
on (pi{Yi). In that direction, we work conditionally on ipi{Yi) = s° = (si, . . . , S2Z+1) 
and denote by Nk{i) = \Yk H Si\ the number of sites of the arc Sj which are occupied 
by Ffc- Write Nk for the family {Nk{i)}i<i<2i+i- Let ii denote the unique index such 
that the singleton Yi_i\Yi C s^. In other words, ii is the unique index i such that 
Ni_i(i) = Ni{i) + 1. Then (pi{Yi_i) results from ipi{Yi) = s by merging the arcs Sj^, Sj^+i 
and Sj;_|_2 (with the convention that indices of arcs are taken modulo 2/ + 1). By iteration, 
we realize that the sequence Nq.j = {Nq, . . . , Ni) determines (Pi{Y)q.j and therefore also 
(p{Y)o:i. Hence it suffices now to check that the conditional distribution of Nq.j given 
Yi only depends on (fii{Yi) — s°, which is straightforward from the dynamics and the 
observation that for every i = 1, . . . , 2Z + 1, the arc Sj has exactly (| s^l + 1)/2 sites which 
are not occupied by Yi. 

We are now able to prove (i). Take t e V^^+i with P {(pi(Yi) = s°, (^(F^+i) = t) > 0. 
Then 

P {ip{Y)o:i = ro,, ifiiYi) = s°, vp(l^+i) = t) 
= J] P iip{Y)oa = ro, \Yi^x, F^+i = y) P {Yi = x, Yl+^ = y) . 

!/e¥>-i(t) 
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Since F is a Markov chain, it follows that for x,y & Cn with P(y; = x, Y^+i = y) > 0, 

P (<^(r)o;/ = ro, \Yi = x, = y)=F (<^(F)o;z = tq, | 1^ = x) . 

Plugging this into the above formula and using the conditional independence of (p{Y)q.i 
and Yi given ipi{Yi), we deduce that for x G (pi'^{s°), 

P((^(F)o, = ro, I Vi{Yi) = s°,<^(r,+i) = t) = F{^{Y)o., = fq, \ Yi = x) . 

Similarly, one sees that the right hand side equals F{(p{Y)q.i = Tq-i \ (pi{Yi) = s°), and 
(i) follows. We turn to (ii) and the description of the transition mechanism. We keep 
the conditioning on ifi{Yi) = s°. Note that Yj+i evolves from Yj by removing uniformly 
at random one of the n — I masses. By identifying Yi with S{Yi), this amounts to 
switching one of the upward steps chosen uniformly at random into a downward step. 
More precisely, under our conditioning, the ith arc Sj is picked with probability 

number of upward steps over (|sj| — l)/2 

total number of upward steps n — I ' 

then one of the upward steps over Sj is selected with uniform probability and changed 
into a downward step. Up to a vertical shift in space, S{Yi) restricted to the arc Sj obeys 
the law of simple random walk conditioned on H_i = |si| (with an obvious modification 
for the last arc S2z+i). Given an upward step over Si is switched, S{Yi-^-i) restricted 
to Si can therefore be regarded as simple random walk conditioned on H^^ = |si|. In 
terms of f{Y), we deduce that (f{YiJ^i) is obtained by first picking the ith arc with 
probability given in then splitting its length according to /i|s.| into three numbers 
ri, r2, corresponding to the first three ladder epochs of simple random walk conditioned 
on if_3 = |sj|, and finally ranking them together with the numbers |sj|, j ^ i, in 
decreasing order. In particular, we realize that for predicting ip{Yi+i) out of (pi{Yi), the 
additional information given by fi{Yi) compared to (p{Yi), namely the location of the 
arcs, is irrelevant. Hence also (ii) holds. □ 

Let us now characterize <f{X). 

Proposition 2.3. f{X) = {f{Xk)}o<k<n is a Markov chain. Its transition mechanism 
from time I < n — 1 to I + 1 is described as follows. 

(a) Conditionally on (p{Xi) = s = (si, . . . , S2(n-0+i) ^ '^2n+i> select an index l out of 
the set of all 3-subsets of {1, . . . , 2(n — /) + 1} according to the law 

P (. ^ k) I ,ix.) ^ s) ^ p,,^,:;,;^:;^^!,),,) . 

(b) Given (p{Xi) = s and l = {i,j, k}, rank the sum r = Si + Sj + Sk together with the 
numbers Sm, £ {1, • • • , 2(r2 — /) + l}\{i,j, k}, in decreasing order to obtain a 
new mass partition. 
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Figure 2: The transition mechanism from ^p{Yi) to ^{¥2), where n = 6. Here, at time 1 
the chain Y is in the configuration state Yi = {1,2,6,8,9}. Then the mass at position 
8 is removed. For the corresponding path, this means that the upward step at time 8 is 
changed into a downward step. The new path 5'(y2) coincides up to time 8 with the old 
path S{Yi) and is then indicated by the dashed fine. The excursion interval [6, 13) is 
broken into three intervals [6,9), [9,12), [12,13). Therefore, v^i(1^2) = (31,82,33,34,35) 



with 31 = [0,1), 32 = [1,6), 33 = [6,9), 34 = [9,12) 35 = [12,13), and ^{Y^) = 
(5,3,3,1,1). 



Proof: From the duality ([2]) it follows that if{X) is obtained by reversing (f{Y) in time. 
In particular, the Markov property carries over from ip{Y) to ip{X). 

It remains to look at the transition mechanism. The step from I = ra— 1 to n is obvious 
from the construction of X and if. Now fix / G {0, . . . ,n — 2}, and let M = 2{n — Z) + 1. 
We work conditionally on (p{Xi) = s = (si, . . . , sm) G "^L-m- By construction, {p{Xi^i) 
is obtained from ip{Xi) by summing up three numbers Si,Sj,Sk, where i,j,k are pairwise 
distinct, and rearranging the sum together with Sm, rn 7^ i,j,k, in decreasing order. 
Write 3° = (si, . . . , 3m) for the partition (fi{Xi), and let u be uniformly distributed on 
{0, . . . , M — 1}, independent of Xi. By the random walk representation and Theorem 
1 of [8J, the law of the cyclically ordered arc lengths (|3i+;^|, . . . , |3jvf-i-jy|) (indices are 
taken modulo M) agrees with the law of the M subsequent ladder epochs of simple 
random walk conditioned on H_m = 2n + l. In particular, the law of (|si+i,|, . . . , |sM-f!/|) 
is invariant under permutation and therefore given by {sa{i), ■ ■ ■ , Sa{M)) , where a is a 
permutation of {1, . . . , M}, chosen uniformly at random and independently of Xi. Note 
that this can also be deduced directly from the fact that Xi is uniformly distributed on 
the space of all /-subsets of {0, ... , 2n}. The probability that Si,Sj,Sk are replaced by their 
sum is given by the probability that the arcs So-i(j)+y, s^-l(^j)_^_^, So-i(;j)_,_^ merge. This is 
the case if and only if the arcs adjoin each other and the singleton Xi^i\Xi is contained 
in that arc which is followed in clockwise order by the other two. More formally, the arcs 
merge if and only if there is a permutation p of the indices i, j and k such that X;+i\X; C 
So-i(p(i))-fi/, and (T~^(p(j)) = a~^{p{i)) + 1, a~^{p{k)) = a~^{p{i)) + 2 (both equalities are 
taken modulo M). Given X;_|_i\X; C So-i(j)_,_j,, the probability that 3o-i(j)_,_y, 3o-i(j)_,_j,, 
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s^-i(^k)+u merge is therefore 

2 1 



M- 1 M-2 
The probabihty that Xi+i\Xi C So-i(j)+,y is 

number of vacant sites in So--i(i)+i/ at time / (sj + l)/2 
total number of vacant sites at time / 2n -\- 1 — I 

Altogether, given ip{Xi) = s, 

P(s,-i(i)+„s,-i(,)+,,s,-i(fc)+, merge) = 2n + 1 - I ) 



X 



M-2' 

which is the probability in (a) in the case / < n — 1. □ 

Theorem 12.11 now easily follows. Indeed, from the last proposition we see that (p{X) 
is equal in law to the skeleton chain {Af^}o<fc<n started from particles of unit mass. 
By the duality relation ([2]), reversing (p{X) in time yields the process ^p{Y), which is the 
state chain of a fragmentation process. 



3 Random binary forest representation 

In this section, we give a second construction of the skeleton chain of the ternary co- 
alescence in terms of random binary forests. The connection between random forests 
and coalescent processes was first observed by Pitman in [18j. In our description, we are 
guided by Chapter 5.2.3 of Bertoin [Tj. 

3.1 Basic definitions on graphs 

We first collect some basic notions on graphs which will useful for our purpose. 

A (undirected) graph is a pair G = {V,E), where V is a finite set and E C {U C 
\^ : |f/| = 2}. The elements of V are called vertices, the elements of E edges. The size 
of a graph is the number of vertices \V\. A subgraph of a graph G = (V, E) is a graph 
H = {V, E') with V dV and E' C E. 

Now let G = {V,E) be a graph. Two vertices v,w are adjacent, if {v,w} G 
E. The degree of a vertex v is the number of vertices adjacent to v. A sequence 
{vi,ei,V2, . . ■ ,Vm,em,Vm+i) such that m > Vi ^ Vj for i ^ j and Ci = {vi,Vi+i} e 

for 1 < 2 < m is called a path, or also a f i-f^+i-path. A cycle is a sequence 
(fi, ei, . . . , Vm, Cm, vi) such that m > 2, {vi,ei,..., Vm-i, em-i,Vm) is a path and e„ = 
{vm,Vi} G E. We say that two vertices v,w are connected, if there exists a v-w-path. 
If there is a v-w-path for any v,w ^ V, we say that the graph G is connected. The 
maximal connected subgraphs of G are its connected components. A connected graph 
without a cycle (as a subgraph) is called a tree. In a tree, a leaf is a vertex of degree 
equals 1, while the vertices of degree greater than 1 are called internal vertices. 
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We are interested in a special family of trees. A binary tree is either a tree consisting 
of a single vertex only, called the root of the tree, or a tree where exactly one vertex has 
degree 2, which we then call the root of the tree, and all the other vertices have degree 
3 or they are leaves. The height of a vertex v in a binary tree is the number of edges of 
the (unique) v-r-path, where r is the root of the tree. If v is not a leaf, then there are 
exactly two vertices w,w' adjacent to v with height strictly bigger than that of v, the 
children of v. We call the pair {{v,w}, {v,w'}} the outgoing edges (from v). Finally, a 
binary forest is a graph such that its connected components are binary trees. The leaves 
or internal vertices of such a forest are then all those of its tree components. 

Observe that a binary forest on N vertices with m tree components has N — m edges, 
{N + m)/2 leaves, and {N — m)/2 internal vertices. 

Remark 3.1. In the literature, a binary tree in our sense is often called a (rooted) full 
labeled binary tree. Here, the term "full" reflects the fact that every vertex other than 
the leaves has two children, and "labeled" stresses that the vertices are distinguishable. 
However, we will use the term "labeled" to indicate a labehng of internal vertices. 

3.2 Dynamics 

Our concern here is to describe the dynamics on the space of binary forests, which will 

lead to another representation of the ternary coalescence. 

As before let iV = 2n + 1. We consider y = {1, 2, . . . , A^} as a set of vertices. Given a 
binary forest on V , we enumerate its tree components according to the increasing order 
of their roots. 

We will assign additional labels to all internal vertices of such a forest. A labeling 
of a binary forest on V with m tree components is a bijective map from the set of 
(A^ — m)/2 internal vertices into {1, . . . , (A^ — m)/2}. A labeled binary forest on V is 
then a binary forest together with a labeling. The set of all labeled binary forests on V 
with m tree components is denoted by J-'(m, A^). Note that T{m, N) is empty if m is 
an even number. 

For every 1 < /c < n, we define a map R : J^{2k — 1, A^) — )■ F{2k + 1, A^) as follows. 
For each r G J^{2k — 1, A^), select the internal vertex with the highest label and delete 
both outgoing edges (and the label, since the vertex is now a leaf). We obtain a labeled 
binary forest with 2A; + 1 trees, which we denote by R{t). 

As the reader may already guess, the map R will be the building block of the fragmen- 
tation mechanism - it breaks the tree with the highest label into three (new) trees. The 
reverse dynamic will correspond the the coagulation mechanism: Out of a binary forest 
with at least three trees, pick one leaf and connect it by adding edges to two distinct 
roots from other tree components. Then, three trees have merged into one (new) tree, 
and the selected leaf has become an internal vertex. Before underlying this procedure 
with randomness, let us analyze the map it! in detail. 

Lemma 3.1. For every 1 < k < n, the map R : J'{2k -1,N) — > T{2k + 1, A^) is 
surjective. More precisely, for every r G J-'{2k + 1, A^), 

|{f G J^{2k -1,N): R{f) ^T}\^{n + k + l)k{2k - 1). 
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Figure 3: From left to right (right to left) one step in the coagulation (fragmentation) 
mechanism is shown. For simplicity, only the labeling of the internal vertices is depicted. 
On the left, the leaf with number 6 on the right is chosen as well as two roots from 
different tree components. They are connected by two edges visualized by the dashed 
lines on the right side. 

Proof: Let r G J-'{2k + 1, A^). In order to construct a generic f G i?^^(r), pick a leaf i 
from r. Write p{i) for the root of the tree component containing i. Then select two roots 
j 7^ j' different from p{i), add the edges {hj'} and label the vertex i with the 

number n — k + 1. Out of three components, we have obtained a new labeled binary tree 
with root p{i), which is part of a forest with 2k — 1 trees. Clearly, this forest is contained 
in i?~^(r). Moreover, different choices of i,j,j' give rise to different forests. To finish 
the proof, note that there are n + k + 1 possible choices for a leaf i, and 2k{2k — l)/2 
possible choices for distinct roots {j,j'}- □ 



Remark 3.2. Applying the map R at most n times destructs a labeled binary forest 
into its single vertices. Due to the recursive structure of trees, this method enables one 
to compute various combinatorial quantities. For example, using |J-'(A^, A^)| = 1 and 
iteratively the identity 

\J^{2k -1,N)\ = {n + k + l)k{2k - 1) \J^{2k + l,N)\ 



provided by Lemma 13. 1[ one obtains for k = 2, . . . ,n + 1 

2^-("+%(2n + l)!(2n-l)!(fc-2)! 
'-^^^^ " " in + k)\ik-l)\i2k-3)\ • 

In the case k = 1, 

|.F(1.7V)|^ ^'"<f>''y^" ^2-»(2n+l)!„!C^. 

[n + Ijl 

where C„ = {2n)\/{{n + l)!n!) is the nth Catalan number. Since there are n! different 
labelings of internal vertices, we deduce that the number of binary trees on V is given 
by 2-"(2r2 + 1)! C„. 
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3.3 From forests to mass partitions 

Denote by the kth concatenation of R, where i?" is the identity map. We randomize 
the input by endowing the space J^(l, N) with the uniform probabihty measure and 
interpret the maps R'' as random variables 

R^:J^{1,N) — ^ J^(2A; + l,iV), k = 0,...,n. 

In words, R^{t) is the forest with 2k +1 tree components which arises from r G J^{1, N) 
by picking the k internal vertices with the highest labels and deleting their outgoing 
edges. By induction, we deduce from Lemma [3.11 that R^ obeys the uniform law on the 
space J-'{2k + 1), for each k. We then consider the random variables 

\R'\^:Hl,N)^Vl_,„ k = 



2n+l' 

2n+l 



where for a tree r G J^{1, N), \R^\^{t) = s = (si, . . . , S2fc+i) G V^n+i the sequence of 
the sizes of the tree components, ranked in decreasing order. 

Turning back to the ternary coalescence, let A*^, k = 0, . . . ,n, denote the skeleton 
chain started from particles of unit mass. Its connection to the sizes of the tree 
components is given by 

Proposition 3.1. The sequence of random variables {|-R"~^|~''}o<fc<n is the state chain 
of the ternary coalescence, that is 

{\R^\\\R--'\\...,\R'\^)^{x;,,...,x:,). 

Proof: For each tree r G J^{1, N), the forest R^{t) has no edges, so \R^\^ = (1, . . . , 1) = 
Xq. Note that given \R^\^ = s = (si, . . . , S21+1) for some 1 < / < n, the mass partition 
|^i-i|4, jg obtained from s by replacing three elements Si,Sj,Sk, where i,j,k are pairwise 
distinct, by their sum. Furthermore, observe that the random variables R^, I < k < 
n, are measurable with respect to the sigma-field generated by RK In particular, by 
Proposition 12. 3[ the claim follows if we show that for every < / < n, for every s = 
(si, . . . , S2(n-i)+i) e 'PL+i and for every 3-subset {i,j, A;} C {1, . . . , 2{n - /) + 1}, 

P = 5^®^-©'= I R^'\ = s) + + ^ 



(2n + 1 - l)2{n - l){2{n - /) - 1) ' 

Take a forest r G J^{2{n — /) + 1,A^). We work conditionally on i?"^' = r. By our 
observation above, is uniformly distributed on the set of (2n + 1 — l){n — l){2{n — 

I) — 1) forests which can be obtained from r in the way described in Lemma IXTl We write 
Ti, • • • , T2{n-i)+i for the tree components of r. For every 3-subset {a, 6, c} C {1, . . . , 2{n — 
I) + 1}, the probability that the leaf i is picked in and the roots are chosen from 
and Tc is therefore 

|r„| + 1 1 
- — ■ X 



{2n + 1 - l){n - l){2{n - I) -1] 
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Hence the probability that i?"^'^^ evolves from r by merging the trees r^, Tb and Tc, 
that is the probability that the leaf i is picked in either r^, Tf, or Tc and connected to the 
roots of the other two components is 

kal + Inl + \Tc\ + 3 



(2n + 1 - l)2{n - l){2{n - /) - 1) ' 

□ 

As a consequence of Theorem 12. H the time-reversed process {|-R'^|~''}o<fc<n is a frag- 
mentation chain with dislocation law fi. 

Remark 3.3. Adapting the proof of Corollary 5.7 in [7] to our situation, we find another 
way to prove Corollary ll.H based on the binary forest representation. Namely, with 
m = 2{n — /) + 1 and s = (si, . . . , Sm) G "^L+i^ there are 

2n + l \ (2n + l)! 



possibilities to partition the set of vertices {1, . . . ,2n + 1} into non-empty disjoint sets 
Ei, i = l,...,m, such that \Ei\ = Si and minii^j < minEj for i < j < m. Without 
labeling internal vertices, the number of binary tree structures which can be attached 
to Ei is Si)\/{{si — l)/2)!. Having chosen a binary tree structure for each E^, there 

are /! possible ways to label the / internal vertices. Recall that the tree components of 
a forest are enumerated in increasing order of their roots. It follows that the number of 
binary forests r G J-'{m, 2n + 1) with tree components such that |rj| = Si is given by 

(2n+l)!/! \ni,s,)\ 



ml 

1=1 



Since i?" ^ is uniformly distributed on J^{m,2n +1), we deduce from Proposition 13.11 
that 

F(r-(s c))- ^^'^ (2n + l)!/! ^ 1^(1,..)! 
^ ' |^(m,2n + l)| ml 

where 7(s) has been defined in Section [T72l Plugging in the values for \J^{m,2n + 1)| 
and \ J^{1, Si)\ from Remark 13.21 results in the expression obtained in Corollary ll.il 



3.4 Encoding forests by paths 

We conclude our discussion of binary forests by illustrating a direct connection to the 
random walk representation described in the preceding section. For this part it is more 
convenient to consider (rooted unlabeled) plane trees and forests. In a plane forest 
vertices are regarded as indistinguishable, but the set of children for each vertex is 
ordered, as well as the set of roots of the different tree components. The ordering induces 
serveral natural enumerations of the vertices. For example, one of them is provided by 
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the order in which the vertices are visited by a depth-first search, see Figure |H More on 
this can be found in Chapter 6.2 of Pitman |19| . 

We will look at (full) binary plane forests. To relate them to the binary forests 
considered above, note that the number of binary plane forests on vertices with k tree 
components is equal to 

W((jV-t)/2)! l^'^-^'l- 

since there are 2'^^"'^^/^ possible orderings of the children of the internal vertices of a 
forest in J^{k,N), k\ orderings of the roots, but neither vertices are labeled nor there 
is an additional identification of internal vertices. Clearly the ternary coalescence with 
a monodisperse initial configuration can also be realized on the space of binary plane 
forests, with the same dynamics. 

There are various ways to code plane trees and forests by discrete functions. For a 
(finite) plane tree 6 on N vertices, one common way is to look at its Lukasiewicz path 
{xi}o<i<N- Denoting hy Vq, . . . , f tv-i the vertices of 6 listed in the order of a depth-first 
search and by k{v) the number of children of vertex v, one defines 

i-i 

= Y,iKv^) - 1), < J < N. 

i=0 

Note that xq = 0, = —1, and 

Xj - xj.i = k{v,_,) - 1, l<j<N. (6) 

It is easy to see that there is a bijection between Lukasiewicz paths and rooted plane 
trees. A sequence of such trees may then by encoded by gluing together the correspond- 
ing Lukasiewicz paths, retaining the relationship ([6]). In other words, the coding of the 
next tree starts if a new minimum is attained. 

Turning to random trees, it follows from Proposition L4 of Le Gall jl5] that a Galton- 
Watson tree with offspring distribution ri{k) = l/2{6o{k) + 62{k)), conditioned to have 
total progeny size A^, is distributed according to a tree chosen uniformly at random 
among the set of all binary plane trees on A^ vertices. Furthermore, the corresponding 
Lukasiewicz path of such a Galton- Watson tree is distributed as the path of simple 
random walk on Z up to time A^, conditioned on = A^ (see Corollary L6 of |15]). 

Putting the observations together, we realize that for an integer < I < n, the 
path of simple random walk up to time A^, conditioned on H_(^2i+i) = N, encodes a 
forest uniformly distributed over all binary plane forests on A^ vertices with 21 + 1 tree 
components. In particular, the sequence of the sizes of the tree components is distributed 
as the sequence of the ladder epochs of the conditioned random walk path, if both are put 
in random uniform order, say. However, we should like to point out that the connection 
is only static in the sense that the dynamics defined on the space of mass partitions does 
not correspond to the dynamics on the space of binary forests as described above. 
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Figure 4: On the left side a binary plane forest on 11 vertices with 3 tree components 
is shown, where the vertices are enumerated by a depth-first search. The corresponding 
Lukasiewicz path is depicted on the right side. The crosses indicate where the coding of 
a new tree starts. 

4 Asymptotics of the ternary coalescence 

Having the concrete realizations at hand, we are now able to investigate asymptotic 
properties of the ternary coalescent process. Let us write X^^^ = t > 0) for 

the coalescence with kernel k, started from the monodisperse configuration (1, . . . , 1) 
consisting of = 2n + 1 atoms of unit mass, and put X'/}^^ = X^^^{Tk), k = 0, . . . ,n. 
The number of particles at time t > is denoted by and the number of jumps 

up to time t by J[^l(t). 

We will consider the space of mass partitions with total mass bounded by 1, 



and the subset Si C Vm of sequences with Yl'i^i ■Si = 1- We equip Vm with the uniform 
distance. The induced topology coincides with that of pointwise convergence and turns 
Vm into a compact space. The /i-distance induces a finer topology. However, if (s„, n G 
N) is some sequence in Vm converging pointwise to s G Si, then the convergence does 
also hold in the /i-sense, as it can be easily deduced from Scheffe's lemma. Therefore, 
on Si all these types of convergence are equivalent. 

We turn to our main result of this section. Recall that the standard additive coa- 
lescent X = (X(t), t G M) is the unique additive coalescent process such that for each 
t G M, X{t) has the law of the ranked sequence ai > a2 > ... of the atoms of a Poisson 
random measure on (0, 00) with intensity measure A{da) = e~^da/ V^na^, conditioned 
on Xli^i 3-1 = 1. We refer to [3] and [12] for background. 

Theorem 4.1. As n ^ 00, the Eii-valued process 



V, 



m 




converges in the sense of finite- dimensional distributions towards the standard additive 
coalescent. 
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Here, the multiplication with is meant element-wise. At first glance the con- 
vergence may look surprising, since the standard additive coalescent has been obtained 
by Evans and Pitman in [12] as a limit of additive coalescent processes as follows. Let 
^[n] _ (^3fN(;(;)^ t > 0) be the stochastic coalescence with binary coagulation kernel 

k{r, s) = r + s, r,s > 0, 

started from the monodisperse configuration with n atoms, each of mass 1/n. Then, as 
n — )■ oo, the time-shifted processes (Xt"'l(t+ (1/2) \nn),t > —(1/2) Inn) converge weakly 
to X. 

However, our convergence result concerns only the finite-dimensional laws. At least 
for the one-dimensional distributions, one might expect a result in this direction if one 
compares the one-dimensional statistics of the skeleton chains of the ternary and the 
additive coalescence X'"''. The states of the latter can be expressed in terms of indepen- 
dent standard Borel variables (see for example (30) in [121), which have a similar tail 
behavior as the hitting time H^. 

Let us briefiy comment on the scaling in the theorem. To obtain a limit for the 
normalized sequence of masses, the number of atoms must be of order ^/n. We refer 
to the forthcoming Lemma 14.21 for a better understanding. As Lemma 14.11 shows, if the 
process A"!^' runs for time t/N^^"^, then the amount of particles has typically reduced from 

to about / 1. Note that when approximating the standard additive coalescent with 
the processes X'"^ starting from n atoms of mass 1/n, the macroscopic picture appears 
at times t+ (1/2) Inn at which there are about i/n/e* particles. Here, roughly speaking, 
the standard Borel law plays the role of the hitting time distribution. Precise statements 
can be found in the books of Pitman [19j, Chapter 10.3, or Bertoin [7J, Chapter 5.3. 

We shall present three different ways to obtain convergence for the rescaled ternary 
coalescence of which we discuss two in details. The first more general method will lead 
to one-dimensional convergence in Proposition 14.11 It relies on the observation that 
the distribution of the hitting time Hk is in the domain of attraction of a stable(l/2) 
law. Then a size-biased reordering is used to construct the limiting mass partition. The 
second method resulting in finite-dimensional convergence (and therefore in the proof 
of the theorem) is more specialized to our situation. It is based on the identification of 
configurations with mass partitions via paths, as described in Section 12.11 Since the two 
methods do not rely on each other, the reader in a hurry may safely skip Section 14.21 
In the closing remark we outline a possible third way to establish finite-dimensional 
convergence, using the random binary forest representation. 

4.1 Number of particles 

In order to relate the behavior of X^^'^ to that of its skeleton chain, we prove a limit 
theorem for the number of particles. As just remarked, it will become clear later why 
we choose the spatial scale factor N~^^'^. 
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Lemma 4.1. For every t > 0, as n ^ oo, 



#[^l(t/iV3/2) 1 

= y - m probabinty. 

N t 



Proof: Using the relation = — 2jt^l(-), the claim will follow once we show 

that 

^ ' - — — ^ in probabihty. (7) 



2 

Remember that J^^\t/N^/^) = max{k e Z+ : N^^'^T^ < t}, where T^. is the fcth coag- 
ulation time given by Tk = X^^Li the rates a{i) = a{i,N) are as in ([1]) and 
61,62,... is a sequence of independent standard exponential variables. Heuristically, 
replacing by its expectation X]i=i'^(0~^! ^^e number of jumps J{t/N^/'^) should 
roughly behave as the maximal k such that N'^^'^'Y^i=i'^i}) '^ — ^- show that 

with the choice kn = n — t^^\/N /2, 



i=l 

where we agree that the sum runs from 1 to the largest integer below /c„. First note that 

^ ^' ^^{N+l-^){N+l-2^){N-2^) 

^""(S (iV-.)(iV-2., )^^(--")- 
Furthermore, some simple computations show that for each e > 0, 

{N -i){N -2iy ~ Jo {N -x){N -2xy^ ^'' ^ 

1 



A^(A^ - 2kr. 



^ +0(n-2+-). 



iV3/2 

Altogether, we obtain ([H]). Moreover, since 

Var [N'/'n^] = N'J^ai{)-' = 0{n-'/'), 

i=l 

we deduce that N^^^Tk„ t in probability. From this ([7]) readily follows. □ 
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4.2 Mass partitions induced by Poisson measures 

We shall now prove one-dimensional convergence of the ternary coalescence process. The 
reader who is not familiar with Poisson random measures and size-biased reorderings is 
invited to consult Section 2.2.3 of Bertoin [7j first. 

Let ^1, ^2, • • • be a sequence of independent copies of H^i. Recall that by Lemma [2?T| 
as / — )• oo, 

P(ei = 2/ + l)~iy5. (9) 

For k G let S2fe+i = + - • -+^2^+1; and denote by 5'(2fc+i.^) ^ random mass partition 
distributed as the rearrangement in decreasing order of ^i/iV, . . . , C,2k+i/N, conditionally 
on S2fc+i = N. As a special case of Corollary 2.2 in [7] we have 

Lemma 4.2. Fix b > 0. Then 5'(2fc+i'^) converges in distribution on Vm as k, n ^ oo 
with k ~ bn^^'^ to the ranked sequence (ai,a2,...) of the atoms of a Poisson random 
measure on (0, oo) with intensity A{da) = brr'^^'^a'^^'^da, conditioned on Xli^i = 1- 

Proof: For k < n, denote by (^i^, . . . ,^2fc+i Af) ^ (2^ + l)-tuple distributed as a size- 
biased reordering of (^i, . . . , ^2fc+i) given T,2k+i = N. It easily follows that for / = 
0, . . . ,n — k, 

P(e,V = 2/ + l) = ^^^^^^^f^^P(6 = 2/ + l|S2,+i = Ar) 

(2fc + l)(2/ + l) _ ^^ P(S2, = Ar-(2/-H)) 

If we fix a G (0, 1), 6 > and let /, k, n tend to infinity with / ~ an, k ~ bn^^^, we obtain 
from ([9]) 

(2fc + lK2/ + l) p ^ 2/ + 1) ~ bn-'/'a-'/%-\ (10) 

Setting Qk = Sir'^k"^, we see again by (|9]) that for — )■ oo, 

(2fc + l)P(ei>^?,.)~l- 

Moreover, since k ~ bn^^'^, we have as A;, n — )■ oo 

Ok 462 
N ^ n ' 

It then follows from the theory of stable laws (see Breiman [10], Chapters 9 and 14) that 



N 



— i- in distribution as k,n ^ oo. 



where ^ = Yl'iZi ^"^^ (^i' ^2, . . .) is the ranked sequence of the atoms of a Poisson ran- 
dom measure on (0, oo) with intensity bn^^^'^a~^^'^da. In particular, ^ is stable(l/2) and 
has a smooth density p{x) = P(<^ G dx)/dx, which is strictly positive on (0,oo). Using 
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/ ~ an, we infer from Gnedenko's local limit theorem (see Gnedenko, Kolmogorov [13J, 
p. 236) that 

P(E2fc = iV-(2/ + l)) p{l - a) 
P(S2fc+i = iV) " p(l) ■ 

Together with (fTOjl this shows 

P (eiV = 2/ + 1) ~ aP^^^bn~'/'a-^/\ 

np[l) 

In particular, C,l^/N converges weakly as n, — )■ oo with /c ~ bn^^^ towards the law 

^p(l-a) ^^_,/2^_3/2^^^ a G (0, 1). 

Now observe that given ^ = 21 + 1 for some / = 0, . . . , n — fc, we have the equality in 
law 

fe,iV5 ■ • • 5 ^2fc+l,Af) = (^l,Af-(2/+l)5 • • • 5 ^2k,N-(2l+l)) ■ 

By iterating the argument from above, we may therefore deduce from Proposition 2.4 
(i) in [7j that the limit law of [^l j^f/N, . . . ,^2k+i n/^) sense of finite-dimensional 

distributions as n. A; — )■ oo with k ~ bn^^"^ is given by the law of a size-biased reordering 
(a^, a.2, . . .) of (ai, a2, . . .), where the latter sequence is as in the statement. Also, we 
have that YllZi — ^ almost surely. From Proposition 2.3 of [7j it then follows that 
the (ranked) random mass partition 5'(2fc+i>^) converges in distribution to the ranked 
sequence of atoms (ai, a2, . . .), conditioned on <j = 1. □ 



For the skeleton chain we derive the following consequence. 



Corollary 4.1. Fix b > 0. If n, k oo with k ~ bn^^^, then (l/N)X^_l converges in 
distribution on Si to the ranked sequence (ai, a2, . . .) of the atoms of a Poisson random 
measure on (0, oo) with intensity A(da) = b7i~^^'^a~^^'^da, conditioned on Yl'iLi ^-i = 1- 

Proof: This follows from Proposition 12.11 together with the last lemma. □ 

Combining the corollary with the weak convergence result for the number of particles, 
we easily obtain one-dimensional convergence. 

Proposition 4.1. Fix t > 0. Then 

l^X^^\t/N'/') 

converges in distribution on E>i to the ranked sequence (ai,a2,...) of the atoms of a 
Poisson random measure on (0, oo) with intensity 

da, a > 0, 



V2 
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conditioned on Yl'ili^i ~ ^- particular, the one- dimensional distributions of the 
process 

t^l;t'[^](e7iV3/2), tGM, 



converge to those of the standard additive coalescent. 

Proof: Let K = n - J^^\t/N^l^). Then X^^\t/N^/^) = X^^l^, so we may show the 
convergence for {l/N)X^^l^. From Lemma [4.11 it follows that as n — )■ oo, 



— ^ — )• —= in probability. 



Furthermore, we know from Corollary 14. 1 1 that if /„ is a deterministic sequence of integers 
with /„ ~ ^Jrdr^ I \f2^ then we have the asserted convergence for {\IN')X'^ . 

It therefore remains to argue that we may replace by the random sequence kn- 
To this end, recall that Vm is a compact metric space, so by Prohorov's theorem (see 
Billingsley |t9j. Section 6) the space of probability measures on Vm is relatively compact, 
and we only have to show convergence on Vm in the sense of finite-dimensional distri- 
butions. Since all involved random mass partitions lie in Si almost surely, this leads to 
convergence in distribution on §i. Denote by xf^^ the ith component of {l/N)X^^l . 
The finite-dimensional convergence follows if we show that for all j E N and Aj > 0, as 
n — )■ oo 



E 



exp - ^ X 



[N] 



i=l 



exp 



[12) 



where (ai, a2, . . .) is distributed as the rearrangement in decreasing order of (a*, z G N) 
under P^, see Proposition 2.4 of [7]. Statement (IT^ is equivalent to 



E 



exp 



X] (AixS^' + . . . + A. 

i=l 



X 



[N] 



E 



exp 



(Aiai + . . . + Aiaj 



■t=i 



and we will show the latter. Denote by / : Vm (0, 1] the function 



(13) 



/(s) = exp - ^ (AiSi + . . . + XiSi) , s = (si, Sa, . . .) G Vm- 



Note that f{{l/N)A:^^^{t)) > f{{l/N)X^^^{s)) almost surely whenever t < s. By ([II]) 
we can find deterministic sequences of integers l~ and such that In^^ti^ ^/nt^^ / \/2 
and the probability of the event {/~ < fc„ < l^} tends to 1 as n — )■ oo. But on this 
event, we have by monotonicity 

The expectations of the outer quantities converge to the right side of f lT^ . This finishes 
the proof. □ 
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4.3 Convergence of ladder epochs 

Aldous and Pitman have shown in [3J that the exponential time change 



= X(-lnt), t>0, 



with -F(O) = (1, 0, . . .) transforms the standard additive coalescent into a fragmentation 
process which is self-similar with index a = 1/2. In [4j, Bertoin has given an explicit 
construction of this fragmentation process in terms of ladder epochs of Brownian excur- 
sion with drift, and our finite-dimensional convergence result for the ternary coalescence 
will be based on this identity. 

Let us introduce some notation. We denote by C[0, 1] the space of continuous real- 
valued paths on [0, 1], endowed with the uniform topology. For an arbitrary path u G 
C[0, 1], its ladder time set is given by 



Since C{u) is a closed set, there exists a unique decomposition of [0,l]\£(a;) into a 
countable union of disjoint (open) intervals. We denote by G{u) the ranked sequence of 
their lengths. By filling up with zeros, we may always interpret G{u) as a mass partition 
in Vm- Note that G{u) G §i if and only if C{uj) has Lebesgue measure zero. 

The construction of the dual fragmentation process F in [4j can be summarized as 
follows. Let e = (e(s), < s < 1) be a positive Brownian excursion. For every t > 0, 
consider the excursion dragged down with drift t, that is et{s) = e{s) — st, < s < 1, 
and its ladder time set C{et), which has almost surely Lebesgue measure zero. It has 
then been shown that the law of (G(et), t > 0) and {F{t), t > 0) coincide. 

In light of our representation of the ternary coalescence in terms of ladder epochs, it 
seems natural to establish convergence of these objects. In this direction, the main step 
is to prove distributional convergence of the underlying random paths towards Brownian 
excursion with drift. 

To begin with, take a process {Jn{t),t> 0) distributed as {J^'^\t/N^/'^), t>0), and 
independently of this a Markov chain {X/}o</<n as defined in Section 12.21 Let us first 
fix t > 0, and write Jn = Jn{t)- Remember that given J„, we may identify Xj^ with 
simple random walk up to time A^, conditioned to end at —{2{n — Jn) + 1), 



where ^i{Xj^) = l{Xj^{i) = 1}. By linear interpolation, we define the corresponding 
continuous random path Sn,t on the unit interval. 




s G [0, 1] : uj{s) = inf u 





SnAs) = 2 E UXjJ + {Ns - [Ns\)^^r,si{XjJ \ - Ns, 0<s<l. 
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We shall now prove convergence of the finite-dimensional laws of the C[0, l]-valued pro- 
cess {N~^^'^Sn,t, i > 0). The limiting object {B^Li, t > 0) is distributed as 

[B^.i, t>0) = {{B'^is) - st-\ < s < 1), t > 0) , (14) 

where B''^ is a standard Brownian bridge on the unit interval. In particular, for each 
fixed t, the distribution of B^Li on C[0, 1] is that of a Brownian bridge from to —t^^. 

Lemma 4.3. The C[0,l]-valued process (^N~^^'^Sn,ti ^ > O) converges in the sense of 
finite- dimensional distributions as n ^ oo to (-B^-i, t > O) . 

Proof: Let us fix t > as above and first prove one-dimensional convergence. For 
< s < 1, define 

Wn{s) = 2 5Z ^^(^") + - [Ns\)iyNs\{Xn) | - Ns, 

y i=o 

y i=o 

We may then rewrite Sn,t as Sn,t = — 

The process Wn{-) is linear interpolation of simple random walk up to time N, con- 
ditioned to end at —1. We deduce from a conditioned version of Donsker's invariance 
principle (see Dwass and Karlin [11]) that {N~^^'^Wn{s), < s < 1) converges weakly in 
C[0, 1] to the standard Brownian bridge B^'". 

Concerning the drift part D^, we let 

[Ns}-1 

so that Dn = 2Dn^ + Dn\ Now fix s G [0,1]. A moment's thought reveals that 
conditioned on J„ = n — A; for some A; G {0, . . . , n}, the random variable Dn\s) follows 
the hypergeometric distribution. More precisely, 

/lNs\\ /N-lNs\\ 

P (D«(s) = j) = ^ ^ /jv^^'^ ' max{0, k + [Ns\ - N} < j < mm{k, [Ns\}. 
As a consequence, under the same conditioning, 

^[Dl!\s)]=k^, Var (/)«(.))< A:. (15) 
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Let kn = n — Jn- Choosing £ > arbitrarily small and n large enough, by the law of 
total probability, 

¥{N-^'^\Dn{s)-2kns\ >e) 
< ^ P (iV-V2|/)(i)(,) _ E[D«(s)| > £/3 I = k) nK = k) 

k=0 

= 0(1), 

where the last line follows from f ITT]) . f lTS]) and Chebyshev's inequality. Since by ([7]), 
N~^/'^2knS converges in probability to t"^s, so does N~^/'^Dn{s). In particular, the finite- 
dimensional laws of (N^^^'^Dni'S), < s < 1) converge to those of (t-^s, < s < 1). 
Moreover, -D„(s) is increasing in s, and a similar computation entails that for A large 
enough, as n — )■ oo, 

P(iV-V2z)„(l)>A)=o(l). 

By Theorem 8.4 of Billingsley [9], we conclude that the distributions of N~^^'^Dn{-) form 
a tight sequence. It follows that {N-^/^Dn{s), < s < 1) converges in probability to 
(r^s, < s < 1). Applying now Theorem 4.4 from together with the continuous 
mapping theorem finishes the proof of the one- dimensional convergence. 

The arguments obviously extend to finite-dimensional distributions. Indeed, the 
bridge term Wn is the same for all t, and the drift term Dn converges in probability, for 
each t. Therefore, finite-dimensional convergence follows again from Theorem 4.4 of [9j. 

□ 



Similar as for discrete paths, we introduce for v G [0, 1] the shift operator 6 on C[0, 1], 
Define H : C[0, 1] — ?■ [0, 1] as the first time when the global minimum is attained. 



uj{s + v) — uj{v) , < s < 1 — f 

uls + v-l) - uj{v) + uj{l) - u{0) , 1-v < s <1 



H{u)) = inf |s G [0, 1] : u{s) = inf w 

Clearly, H is not continuous on the whole space, but it is so restricted to the subset 
of paths that uniquely attain their minimum. It is well-known and also implied by 
the forthcoming Lemma 14.41 that the distribution of B'^Li is fully supported on this 
subset. Further, the shift operator is continuous as a map 6 : C[0, 1] x [0, 1] — )■ C[0, 1], 
9{u,v) = 9yU. Setting 9hi^ = Oh{uj)UJ, it then follows from the above lemma and the 
continuous mapping theorem that for n — )■ oo, 

{N-'/'eHSn,t, t > 0) ^ {9hB^:^, t > 0) 
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in the sense of finite-dimensional distributions. Recall < s < f ) = {B^^{s) — 

st~^, < s < 1), where B^^ is a Brownian bridge (the same for all t). Denoting by e a 
standard Brownian excursion, it has been proven by Vervaat in [20J that 

OhB'"^ = e. 

Since 6u o 0v = for w = u + f [mod 1], we have Oh = Oh ° Oy pathwise for every 
Q < V < 1. Therefore, if /i denotes the almost surely unique instant when attains 
its minimum, 

OhB\':, = OhoO^ [B^' - st-\ < s < 1) 
= Oh {OhB'"' - sr\ < s < 1) 

= Onet-^. (16) 

Here, as above, et-i{s) = e{s) — st~^ is the Brownian excursion dragged down with drift 
t~^. Since e^-i attains its minimal value almost surely at the endpoint, we have proven 
the following 

Corollary 4.2. In the notation above, {^N^^^'^OnSn^t^ t > O) converges in the sense of 
finite- dimensional distributions as n ^ oo to (e^-i, t > 0). 

The convergence of the ternary coalescence is now easy to establish. As last prepara- 
tion, let us recall a technical result. Call a point x G [0, 1] a local minimum of a; G C[0, 1], 
if there exists 5 > such that for all y G [max{a; — 5, 0}, min{a; + 5, 1}], uj{x) < uj{y). 
The following statement is true for all real t. 

Lemma 4.4. With probability one, all local minima of {£t{s), < s < 1) are distinct. 

Proof: By (fT6|) . we may show the statement for {B^'^{s), < s < 1) instead. Since for 
the time-reversed process, it holds that 

{B^{1 - s) - (1 - s)t, < s < 1) = (B'^is) + st-t,0<s<l), 

it suffices to show that for some 1/2 < r < 1, (B'l^^s), < s < r) has almost surely 
distinct local minima. However, if J> denotes the filtration generated by the canonical 
process x. on C[0, 1] up to time r < 1, Q the law of B'l'^ and, for a moment, P Wiener 
measure and p the Gaussian transition kernel, it is well-known that Q is locally absolute 
continuous with respect to P, 

Pi{0,-t) 



Since the local minima of Brownian motion on [0, 1] are distinct almost surely (see for 
example Theorem 2.11 in the book by Morters and Peres [IS]), the lemma is proven. □ 
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Proof of Theorem I4.lt In view of Bertoin's result in [4j, the claim follows if we show 
that the finite-dimensional laws of 

t^lxWi^t/N^/^), t>0, 

converge to those of {G{et-i), t > 0). Remember the map if constructed in Section [2?T] 
sending configurations to mass partitions. With J„(t) = J^^^t/N^^^) defined as above, 
we have already seen that 

Let t > 0, and assume that conditionally on 

j;^V{^Mt)) = (Sl, • • • , S2(n-J„{t))+l), 

where Nsi G {1, 3, 5, ... , A^} with ^ Sj = 1. Then by construction of both (f, G and the 
linear interpolation, 

G{N^^^'^6HSn,t) = {gi, ■ ■ ■ ,g2{n-J„{t))+l), 

with (yfj + = s j for all i. Thus, we are left with showing finite-dimensional convergence 
of {G{N-^/^eHSn,t), t > 0) to (G'(et-i), t > 0). It is easy to check that G : C[0, 1] ^ Vm 
is continuous on the subset of those paths which attain their local minima at unique 
points. By Lemma [4.4^ the distribution of e^-i assigns mass one to this subset. There- 
fore, Corollary 14.21 and the continuous mapping theorem yield convergence of the finite- 
dimensional distributions on Vm, and since G{et-^) G Si with probability one, we also 
obtain finite-dimensional convergence on Si. □ 

Final remark 

For the proof of theorem we made use of the random walk representation. Let us 
point out another possibility to derive convergence, using the random binary forest 
representation. Following the construction in Section [21 the state chain of the ternary 
coalescence starting from particles of unit mass can be obtained in reversed time 
by deleting successively pairs of outgoing edges of a random tree uniformly distributed 
over all binary plane trees on vertices. Such a random tree can be seen as a Galton- 
Watson tree with offspring distribution ^{k) = \{So{k)+62{k)), conditioned to have total 
population size A^. One finds oneself in the setting of Theorem 23 (in the sublattice case) 
of Aldous P]. In particular, if r^^l denotes the uniform binary plane tree on A^ vertices, 
where mass 1/A^ is assigned to each vertex and the edges are rescaled to have length 
1/\/N, then r^^' converges weakly as A^ — )■ cxd to the Brownian continuum random tree 
(CRT) introduced in |JJ. The deletion of edges in the binary plane tree corresponds to 
cutting the skeleton of the CRT into subtrees. It has been shown by Aldous and Pitman 
in [3] that if for each — oo < t < oo, the skeleton is split at marks set by a Poisson process 
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with rate e~* per unit length, then the process of the ranked masses of these subtrees 
evolves as the standard additive coalescent. It should be possible to approximate the 
Poisson process of marks on the CRT by the process of deleting edges from the binary 
plane tree. This would then lead to another proof of Theorem 14.11 
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